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Abstract 

In  the  analysis  of  a  packetized  communication  system  such  as  a  slotted 
AIDHA,  one  needs  to  deal  with  discrete-time  point  processes  and  related 
queues.  We  examine  how  Poisson  processes,  M/M/1  queue,  M/G/«>,  etc.,  should 
be  translated  in  their  discrete-time  versions.  We  then  give  a  new  definition 
of  processor-sharing  and  its  properties.  Correlation  functions  of  input  and 
output  processes  in  an  M/G/°°  are  also  studied. 
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I .  INTRODUCTION 


In  the  analysis  of  a  packet ized  communication  system,  one  needs  to 
deal  with  discrete-time  point  processes.  One  may  quickly  conclude  that 
most  results  and  formulas  known  for  continuous -time  processes  can  be  trans¬ 
lated  into  their  discrete  time  versions  in  obvious  manner.  Unfortunately, 
this  turns  out  not  to  be  the  case.  The  main  problem  lies  in  that  a  Poisson 
process  lias  two  discrete  time  analogs,  one  is  a  Bernoulli  sequence,  and  the 
other  is  a  Poisson  sequence.  A  discrete  time  version  of  an  M/M/1  system 
holds  for  a  Bernoulli  arrival  sequence,  whereas  an  lVG/°°  can  be  translated 
only  for  a  Poisson  sequence  as  will  be  discussed  in  Sections  II  and  III. 
Important  results  obtained  for  an  M/G/l  system  under  processor  sharing  (PS) 
cannot  be  interpreted  for  the  round-robin  (RR)  scheduling,  despite  the  fact 
the  notion  of  the  PS  scheduling  was  originally  derived  from  RR.  Thus  we 
introduce  a  new  notion  of  processor  sharing. 

This  note  is  therefore  intended  to  examine  some  fundamentals  issues 
of  discrete  time  point  processes  and  related  queueing  systems.  Applications 
of  these  results  to  modeling  and  performance  analysis  of  a  slotted  ALOHA 
SYSTEM  will  be  discussed  in  subsequent  reports . 


II.  A  POISSON  SEQUENCE  WITH  RANDOM  DELAY:  SYSTEM  M/G/°° 

Definition  2.1:  A  discrete -time  point  process  Xk  is  called  a  Poisson 
sequence  if  {Xk:k=0,l,2, . . . }  are  i.i.d.  with  Poisson 
distribution  of  constant  mean  X. 

Consider  a  point  process  X(t),  and  let  us  define  Xk  as  the  number 
of  points  (arrivals)  in  the  interval  (k-l)T<t<kT.  Clearly,  if  X(t)  is 
a  Poisson  process  with  rate  m,  then  Xk  is  a  Poisson  sequence  with  mean 
X=pT.  For  example,  the  number  of  messages  arriving  into  a  slotted  ALOHA 
channel  from  a  large  number  of  terminals  can  be  accurately  characterized 
by  a  Poisson  sequence. 


Random  Delay 

Suppose  that  each  element  (e.g.,  message  or  job)  of  a  discrete-time 

point  process  is  passed  through  a  "random  delay"  (RD)  as  depicted  in  Figure 

2.1.  Each  element  (message)  stays  in  one  of  the  infinite  boxes  for  a 

random  amount  of  time  a.  The  delays  of  the  individual  messages  are  chosen 

00 

independently,  but  from  a  corunon  distribution  (f(d)},  l  f(d)=l.  Backlogged 

d=0 

messages  in  a  slotted  ALOHA  channel  with  random  access  scheme,  for  example, 

can  be  modelled  in  this  manner. 

Random  Delay 


input 

sequence 
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Figure  2.1 


The  notion  of  random  delay  is  equivalent  to  a  service  station  which 
has  infinite  servers  in  parallel.  In  this  case  f(d)  corresponds  to 
service  time  distribution.  The  notion  of  infinite  servers  is  applicable 
even  to  a  finite  set  of  parallel  servers,  provided  the  number  is  sufficient 
that  no  queue  is  formed.  For  example,  a  collection  of  user  terminals  in 
a  finite  population  model  of  an  ALOHA  SYSTEM  can  be  characterized  by  RD, 
in  which  f(d)  is  the  probability  of  the  terminal  think  time  d. 


For  a  Poisson  arrival  sequence  with  rate  A,  the  probability 
that  there  are  n  jobs  (messages)  in  the  RD  at  time  k  is 


p(n;k)  JJZL 


k-1  r  )n 
X  l  FC(d)> 


(  k-!  r 
exp  <-X  l  FL(d) 
d=0 


(2.1) 


c 

where  F  (d)  is  the  complementary  distribution  of  delay 
d: 

F^(d)  =  Prob[delay>d] 


“  I  f (i)  (2.2) 

i>d 

In  the  limit  k-*»,  we  have  the  following  equilibrium 
distribution 

n 

p(n)  =  e"P  (2.3) 

where 


p  =  AE[d] 


(2.4) 


Furthermore,  the  ^ut  (departing)  process  from  RD  is  also  a  Poisson 
sequence  with  mean  X . 

Proof.  See  Appendix  A. 

Different  Classes  of  Messages  or  Terminals 

Suppose  that  there  are  R  different  classes  of  messages  (or  terminals) 
with  different  delay  (or  think  time)  distributions  Fr(d)=  £  fr(i)  r=l,2,.. 
Let  us  denote  the  system  state  by  vector 

2.  ~  »^2  >  •  •  •  >n^] 

where  n.  is  the  number  of  jobs  in  RD  at  a  given  time.  If  class-r  messages 
arrivals  are  independent  Poisson  sequences  with  rates  Ar,  the  probability 
of  the  system  state  at  time  k  is 


P (n ; k)  =  n  p  (n  ;k) 

r=l  r  r 


(2.5) 


where  p  (n  ;k)  takes  the  form  of  (2.1).  In  the  limit  k-*» 


Rp’V 

PCs)  =  n  -J-r  e_Pr 
j=l  V 


(2.6) 


The  corresponding  probability  generating  function  (p.g.f.)  is  given  by 


,  ,  r  R  nJ 

(z)  =  E  Hz*  = 

Lr=l  r  J 


(  R 

exp  <  -p+  l  p  z 
r-1  r  r 


(2.7) 


where 


R 

P-  l 


ilpr  =  J/W 


(2.8) 


'tv*  v 


Thus,  the  p.g.f.  of  the  marginal  distribution  of  the  total  number  of  messages 
in  RD  is 

Q(z)  =  exp{-  p+pz)  (2.9) 

k 

which  yields  the  distribution  of  n  *  J  rij.  : 

r=l 
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Definition  3.1:  A  discrete  time  0-1  valued  process  is  called  a  Bernoulli 
sequence  with  parameter  X  if  (Xk;k=0,l,2, . . . }  are  i.i.d. 
with  distribution  P[Xk=l]=X  ,  and  P[Xk=0]=l-\,  for  all  k. 

Let  T\  be  the  number  of  zero  runs  (which  hereafter  we  call  "run-length") 
between  the  (i-l)bt  one  and  the  1th  one  in  {Xk>.  Then  (T^  are  i.i.d.  with 
geometric  distribution. 


PfT^t]  =  A  (l-X)1 
Thus  the  average  run- length  is 


£-0,1,2,...  (3.1) 


(3.2) 


For  example,  the  sequence  of  messages  from  a  user  terminal  is  often 
modelled  as  a  Bernoulli  sequence,  and  the  variable  T  is  often  called  the 
"thinking  time".  In  a  discrete-time  queueing  system,  the  Bernoulli  process 
plays  a  role  equivalent  to  that  of  a  Poisson  process  in  a  continuous -time 
queueing,  system. 

Definition  3.2:  A  discrete-time  queueing  system  is  denoted  by  M/M/1  if 

its  message  (job)  arrival  sequence  is  a  Bernoulli  sequence, 
the  message  lengths  (job  service  times)  are  i.i.d.  with  a 
geometric  distribution,  and  there  is  only  one  channel 


Definition  3.3:  A  queueing  (scheduling)  discipline  is  said  to  be  a  work- 

conserving  discipline  [8]  if  (1)  the  lengths  of  each  message 
are  not  affected  by  the  queue  discipline;  (2)  the  queue 
discipline  cannot  take  advantage  of  possible  knowledge 
about  message  arrival  times  or  message  lengths;  and  (3) 
the  channel  should  not  be  idle  if  there  are  some  messages 
to  be  sent. 

Because  of  the  so-called  memoryless  property  of  geometric  distribution, 
the  state  of  the  system  M/M/1  under  a  work -conserving  discipline  can  be 
defined  simply  by  n,  the  number  of  messages  in  the  system  (either  being 
transmitted  over  the  channel  or  waiting  in  the  queue) .  Suppose  that  the 
message  generating  probability  is  dependent  on  the  system  state,  and  we 
denote  it  by  X  (n)  [message/slot  time].  Similarly,  assume  that  service  rate 
(or  channel  speed)  is  dependent  on  the  system  state  and  denoted  by  C(n) 
[packets/ slot  time] .  If  the  number  of  packet  length  L  in  a  message  is 
geometrically  distributed  with  parameter  r  (i.e.,  the  average  message  has 

__  1  _v* 

the  length  L=  -p—  [packets]) 

P[L=£]  =  r(l-r)£  £=0,1,2,...  (3.3) 

we  obtain  the  following  result  by  solving  a  discrete  time  version  of  the 
birth-and-death  process: 

Property  3.1:  The  equilibrium  state  distribution  of  n  in  the  M/M/1  under  a 
work-conserving  discipline,  if  it  exists,  is  given  by 

p(n)  =  pWAfnjDtn)!?1 


(3.4) 


A(n)  =  n  X(i-l) 
i=l 


(3.5) 


n 


D(n)  =  n 


i=i^T 


(3.6) 


Note  the  existence  of  the  equilibrium  distribution  (i.e.,  the  condition 
for  system  stability)  is  equivalent  to  the  condition  that  the  p.g.f . 


Q(z)  =  P(0)  l  A(n)D(n)(Fz)n 
n=0 


(3.7) 


is  analytic  inside  the  unit  disc  |z|=l. 

The  unknown  constant  is  given  by  setting  Q(l)=l: 


P(O)'1  =  l  AfoPOOE11 
n=0 


(3.8) 


From  the  above  result  we  can  also  show 


J 


J 


Property  3.2:  The  output  (departure)  process  with  time  "reversed"  is  a 
Bernoulli  process  with  rate  X(n),  where  n  is  the  "current" 
state  of  the  system.  Thus  if  the  arrival  is  a  homogeneous 
Bernoulli  process  with  rate  X,  so  is  the  output  sequence. 

An  exact  analogy  between  continuous -time  and  discrete-time  queueing 
systems  breaks  down,  however,  because  Bernoulli  processes  do  not  possess  the 
reproducing"  property  under  addition.  Consider  a  time  interval  of  arbitrary 
length  [time  slots)  and  consider  m  independent  Bernoulli  sequences.  The 
number  of  arrivals  n.  from  the  kth  stream  during  this  interval  is  the  binomial 
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distribution  with  mean 


p(nk]  - 


Nnk 


(3.9) 


which  has  the  p.g.f. 


^(z)  «  {l-Xk(l-z)} 


N 


(3.10 


Thus  the  total  number  of  arrivals  fran  all  to  m  sources  has  the  p.g.f. 


m 


m 


G(m)(z)  =  n  G,  (z)  =  n  (1-X,  (l-z)}N 
k=l  K  k=l  k 


(3.11. 


If  in  particular  Xk  =  ^  for  all  k=l,2,...,m 

G^\z) = {i_  i  (i.z)r 


(3.12] 


In  the  limit 


lim  G^(z)  =  e"X(1"z)N 
m-*» 


(3.13) 


Thus  the  superposition  of  many  independent  Bernoulli  sequences  with  saj 
rates  converges  in  the  limit  to  a  Poisson  sequence. 


0- 


IL _ PROCESSOR-SHARING  and  general  message  length  distribution 

Definition  4.1:  We  define  processor -sharing  in  a  discrete  time  queueirg 
system  as  a  dispatching  algorithm  which  selects,  in  every 
slot  time,  one  of  n  messages  in  the  system  randomly  and 
equally  likely. 

For  example,  we  can  approximate  a  slotted  ALOHA  channel  I  y  a  processor 
shared  server  with  queue -dependent  rate  C(n).  (How  to  determine  the 
function  C(n)  for  a  slotted  ALOHA  channel  will  be  discussed  in  details  in 
a  forthcoming  report.)  Note  that  the  notion  «f  processor -sharing  in  a 
discrete-time  system  is  different  from  the  so-callod  round-robin  (RR)  sched- 
uling.  When  the  time  quantum  or  slot  time  approaches  zero,  the  notions  of 
RR  and  PS  defined  above  are  both  reduced  to  the  familiar  notion  of  PS  in  con¬ 
tinuous-time  queueing  theory.  We  can  establish  the  following  important  result 

on  PS.  The  corresponding  result  in  continuous-time  queueing  system  is  discussed 
in  many  recent  articles  (e.g.,  [1,2]). 

Property  4.1;  Let  a  message  sequence  be  characterized  by  a  Bernoulli  sequence 
with  rate  A(n),  and  let  the  distribution  of  the  number  of 
packets  in  a  message  be  a  general  distribution  FT  (£)  with 
mean  L  [packets] ,  i .  e . , 


L  =  l  Fh*) 

£=0  L 


(4.1) 


If  the  queue  discipline  is  PS  and  processing  rate  is  C(n)  [packets/slot  time] 
then  the  equilibrium  distribution  of  n  is  given  by  the  formula  (3.4): 


P(n)  =  p(0)A(n)D(n)En 


Proof :  Let  $(z)  be  the  p.g.f,  or  z-transform  of  the  message  length 
distribution; 


$(Z)  =  j0zl(FLCi)-FL(i-1)}  (4.3) 

which  is  a  discrete  analog  of  the  Laplace -Stieltjes  transform. 
If  <J>(z)  can  be  written  as  a  rational  function  of  z: 


(4.4) 


we  can  obtain  the  following  expansion 


where  q  is  the  degree  of  the  polynomial  Q(z),  and  {ct^1}  are  the  character¬ 
istic  roots  of  Q(z)=0.  The  coefficients  (a^  and  {b^  satisfy 

ai  +  =  1  0<i<_q-l  (4.6) 


and 


bq  "  1  (4.7) 

The  representation  (4.5)  is  schematically  shown  in  Figure  4.1,  which  can  be 
viewed  as  cascaded  geometric  distributions*  with  parameters  c1,a2,...,o  . 

This  expansion  is  a  discrete  analog  of  Cox's  representation  [3]  of  a  general 
service  time  distribution  whose  Laplace -Stieltjes  transform  is  a  rational 
function  of  the  Laplacian  variable  s .  Appendix  B  discusses  illustrative 
examples  of  the  expansion  (4.5). 

We  define  the  system  state  by  the  vector  n={n^,l<j<q}  where  n^  is  the 

number  of  messages  in  the  system  which  are  currently  in  the  jth  (fictitious) 

geometric  server  in  the  representation  of  Figure  4.1.  We  define  n(j")  as  the 

state  obtained  by  reducing  its  j^1  component  n.  by  one,  and  n(j+)  is  obtained 

th  J 

u/  adding  unity  to  the  j  component.  Similarly  the  transition  (state  n(j  +  ,i" 
state  takes  place  if  one  of  the  elements  moves  to  the  i^1  component. 
Let  P[n; k]  be  the  probability  that  the  system  is  in  state  n  at  slot  time  k. 

Note  that  this  distribution  takes  a  slightly  different  form  from  (3.3): 

PM)  =  ai"1(1-ai) 


£=1,2,3, . . . 


To  calculate  P[n;k+1] ,  we  note  that  the  system  equation  can  be  written  as  follows 


P[n;k+1]  =  P[n;k] (l-X(n) -  I  ]ma.C(n)}  +  P[nU');k]A(n-l) 

+  n^l  ^  ptH(J+);k]  (n.+l)a.b.C(n+l) 
j=l  J  J  J 

*  n  (4.8) 


where 


a.  =1  -  a. 
J  J 


(4.9) 


The  steady  state  distribution,  if  it  ever  exists,  must  satisfy  the  balance 
equation  (4.8)  in  which  we  set  Ji55{n;k+l]=J™P[n;k]=P[n] .  Then  after  a 
rather  extensive  manipulation  we  can  show  that  the  following  recurrence 
relation  is  a  sufficient  condition  for  the  global  balance  equation. 


n. 

iraiC(n)P[n]  «  e.A(n-l)P[n(k')]  (4.10) 


for  all  i=l,2,...q  and  for  all  states  n.  The  equation  (4.10)  is  a  discrete 
system  equivalent  of  the  local  balance  equation  or  individual  balance  equation 
discussed  in  [4, 5, 1,2].  In  (4.10)  the  new  quantity  e^  is  the  probability  that 
a  message  ever  reaches  the  kth  (fictitious)  service  stage: 

k-1 

ek  =  .^ai  (4.11) 

By  using  the  recurrence  equation  repeatedly  we  obtain 

P[n]  =  P(0)A(n)D(n)n!  n  (4.12) 

i=l  ni'  \«i  / 

Then  the  marginal  distribution  of  n=n1+n2+. .  .+n  is  given  by 

p(n)  =  p(0)A(n)D(n)Ln  14.131 


The  outpv  (departure)  process  from  the  PS  system  has  the  same  property 
as  M/M/1  which  was  stated  as  Property  3,3. 


Different  Classes  of  Messages 

The  above  result  can  be  generalized  to  messages  with  different  classes. 


Property  4.2:  Let  R  independent  message  streams  be  characterized  by 
Bernoulli  sequences  with  rates  f  X(n),  where  f  is  the 

r\  1  T* 


th  K 

fraction  of  the  r  class,  £  f  =1 .  Let  the  message  length 

r-1  r 

distributions  be  Fr(4),  r=l,2,...,R.  If  the  channel 
(or  server)  adopts  the  PS  discipline  with  rate  C(n) ,  the 
joint  distribution,  P t f } ] ,  in  equilibrium  is  given  by 

P[{nr>]  =  P(0)A(n)D(n)n!  II  jjiy-  (Cfp  r  (4.15) 


where  L^  [packets]  is  the  average  length  of  the  class -r  mes¬ 
sage: 

00 

L  =  l  (1-F  (£)}  (4.16) 

r  A=0  r 

Hence,  the  distribution  of  the  total  number  of  messages  in  the 
system  is 


P(n)  =  P(0)A(n)D(n)Ln 


(4.17) 


V.  CORRELATION  FUNCTIONS 


In  this  section  we  discuss  some  statistical  properties  of  discrete  time 
processes  observed  in  a  system  like  slotted  ALOHA.  One  of  the  crucial 
assumptions  usually  made  in  the  analysis  of  such  system  is  the  Poisson 
sequence  approximation  of  the  total  traffic  into  a  channel.  As  the  first 
step  toward  a  more  accurate  characterization,  correlation  functions  of 
the  related  processes  are  examined. 

Let  {X^}  be  a  stationary  discrete  time  point  process  and  be  passed 
into  the  RD  (random  delay  device)  and  we  denote  the  output  sequence  by 

<y.  ^e  define  autocorrelation  and  cross -correlation  functions  as  in  the 
usual  way: 


Rxx(M)  A  e[ (x^-X)  (x^-X) ] 

(5.1) 

Rxy(k,£)  A  E( (X^-X) (Y^-Y) ] 

(5.2) 

etc. 


Property  5.1:  For  a  wide -sense  stationary  sequence  {Xk},  we  have 


R^CM)  = 


(5.3) 


and 

ryy{M)  = 


where  {fk}  is  the  probability  that  the  delay  length  is  k 
[time  slots],  and  «k  is  Kronecker’s  delta. 

Proof:  See  Appendix  C. 


Equation  (5.3)  has  a  resemblance  to  the  well  known  formula  for  the 
cross -correlation  between  a  linear  filter  input  and  output,  in  which  the 
impulse  response  of  the  filter  is  {fk>.  Equation  (5.4)  is  different,  how¬ 
ever,  from  the  corresponding  expression  for  the  linear  filter  output,  which 
does  not  possess  the  second  term  of  (5.4). 

For  example,  if  the  input  is  a  Poisson  sequence  with  rate  then 
Rxx(k,A*)»A6k£,  •  Therefore, 


RyyCM)  -  Af£_k 


(5.5) 


and 


RyyCk.l)  »  X  l  fk.k.fjt.k,»A«k,l 


"  CS.6) 

which  is  consistent  with  the  earlier  result  that  {YkJ  is  also  a  Poisson 


sequence . 


As  the  second  example,  consider  the  case  where  is  a  Bernoulli  sequence 
with  rate  X.  Then  X=X  and 

RXX(k’A)  =  X(1'^6k,iL  (5-7) 

Thus 

RyyCk^)  =  -X  £tfk-lf.£Jl_k.+*<5k>£  (5.8) 

which  is  sufficient  to  conclude  that  is  neither  Bernoulli  nor  Poisson. 
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APPENDIX  A 


Proof  of  Property  2.1 


The  following  proof  is  essentially  a  discrete  analog  of  the  result  known 
for  M/G/°°  [7] .  Consider  a  message  which  has  arrived  by  slot  time  k.  The 
uniformity  of  Poisson  process  indicates  that  the  arrival  time  j  is  uniformly 
distributed  over  k  slots.  The  message  will  be  still  in  RD  at  slot  time  k 
(-J)  with  prob  P[delay>_k-j+l]=l-F(k-j)=F^(k-j) .  The  probability  that  a 
message  is  still  in  Id)  at  slot  k  given  the  condition  that  the  job  lias  arrived 
by  that  slot  is  therefore 

P  “  £  l  FC(k-j)  =  £k£VCd)  (A-l) 

J  =1  d=0 

Clearly 


lim  p 
k-*» 


(A- 2) 


Suppose  that  n  messages  have  arrived  by  slot  time  k.  Since  their 
delay  times  in  the  RD  are  i.i.d.,  the  probability  that  i  messages  are  still 
in  RD  at  time  k  is  given  by 

(i)  P1(l-p)n"1  0<i<n  (A- 3) 

For  a  Poisson  arrival  sequence  with  rate  X ,  the  number  n  is  Poisson 
distributed  with  mean  Xk: 

(MQn  (A-4 J 

n!  ' 

Thus  the  probability  that  there  are  i  messages  in  RD  at  time  k  is: 


By  substituting  (A-l)  into  (A-5),  we  obtain  (2.1). 

With  p  defined  by  (A-l),  1-p  represents  the  probability  that  an 
arbitrary  message  which  has  arrived  by  slot  k  has  departed  already. 
Therefore,  the  total  number  of  messages  j  that  have  departed  by  the  end  of 
slot  k  is  obtained  by  substituting  1-p  for  p  in  (A-5) : 

ttj  CkJ  -  AkQjP)}j  e~>k(1-P)  (A-6) 

k-1 

which  is  again  a  Poisson  distribution  but  with  mean  A  J  F(d).  We  can  show 

d=0 

that  the  number  leaving  RD  at  the  end  of  (k-l)st  slot  is  independent  of  the 
number  of  messages  left  at  slot  time  k.  This  independence  of  the  past 
departure  process  and  current  state  of  the  system  is  important.  It  follows 
that  the  number  of  departures  in  different  slots  are  statistically  independ 
ent.  Thus,  we  can  show  that  the  number  of  departure  in  slot  k  is  Poisson 
distributed  with  mean  XF(k-l).  In  the  limit  J;im,\F(k-l)=X . 


APPENDIX  B 


On  the  Expansion  of  (4.cj 

The  geometric  stage  representation  of  a  general  discrete  distribution 
involves  the  formal  use  of  complex  transition  probabilities,  since  the 
characteristic  roots  ou1  of  $(z)  can  take  on  complex  values.  Of  course 
we  can  choose  a  different  type  of  expansion  in  which  all  parameters  are 
real.  This  involves  in  general  serial -parallel  combination  of  geometric 
distributions  and  will  not  be  given  in  such  compact  form  as  the  expansion  (4.5) 
Notp  also  that  we  do  not  really  need  to  find  expansion  parameters  {a^}  {a.} 
and  {bi>.  The  final  result  (4.13)  does  not  include  these  parameters.  In 
this  sense  the  representation  (4.5)  is  analogous  tc  the  Fourier  series 
expansion  or  general  orthogonal  expansion  used  in  signa  analysis. 

To  give  the  reader  some  feel  of  the  representation  (4.5),  we  will  show 
illustrative  examples  in  which  expansion  coefficients  are  easily  found. 

(1)  The  distribution 

f(k)  =  (k^iX)  (i-a) qak  k*0,l,2, . . .  (B-l) 

is  a  variant  of  the  negative  binomial  or  Pascal  distribution.  Its  p.g.f. 
is 

t(2)  =  Jof(k)zk  ■  (ra)  “  a«) 

Thus  the  expansion  (4.5)  is  simply  obtained  by  setting  all  a^'s  to  be 
unity  and  a^=0(2=. .  .=a^=a . 

(2)  The  Poisson  distribution 


has  the  corresponding  p.g.f. 


$(z)  =  eX(z-1)  (B-4) 
This  is  clearly  not  a  rational  function  of  z.  If  we  allow,  however, 
the  formal  passages  to  the  limit  we  can  represent  $(z)  as 


$(z)  = 


lim 

q-K» 


(B-5) 


See  [6]  for  an  example  in  which  probability  parameter  ou 's  are  complex 


j 


APPENDIX  C 


Derivation  of  Property  5.1 


We  write  the  output  process  as 


Y£  =  l  2f0) 
j=0  £ 


(C-1) 


where  z.(j)  is  a  contribution  fran  X.  to  the  value  Y„,  i.e.. 
x,  J  Z’  * 


Then 


Therefore 


2^0)  =  n  with  ^  °lnlxj  (C-2) 


E[XkY,]  -  .[oF.[XkZ,(j)] 


X, 


J  J  /XA  n  X. -n 

"  j=0  n-0riXk'n  /f^-JCl'fH)  P[Xk’Xj- 


X,-l 


2  p[xk’V 


.i0{RXxW)+X 


(C-3) 


R^y(k,£)  =  E[XkY^J  -  XY  =  .I^^)Q(Ck>j)Fj2,_j  (C-4) 


which  proves  (5.4). 


From  (C-l)  we  write  for  k^£, 


E<W  = E 


I  2k(i)  I  zoCO 
i=0  K  j=0  *  J 


UL 


l  £fk-iEtXiVf*-i  +  H  l  m  E 

1 J  J  J  imn 


X. !  X . -m-n 

1 _ rin  fn  n  r  r  X  1 

m!n!  (X.  -m-n)  !  rk-i  Z-i1-1"1^-!  £-iJ 


*  I  pk-iWA-r 

(itf) 


CC-5) 


Hence  for  k/Jt 


RVY(k’t)  =  I  pk-AxCi'i^-r1  pk-A-i 


CC-6) 


For  k=£,  (C-5)  should  be  modified  as 

2t  _  r  rr  n  r\r  -*r  n  /•  r»  r»  2T 


EtYkJ  =  l  lf k-iEWfk-i  +  X  E 

i  J  J  J  im 

W) 


XA  X.-m 

)  S-i^k-i*  1 


(C- 


Then  writing  m  =  m(m-l)+m,  we  have 


-  l  l£k-i+X  £fk-i 

i  J  J  i  i 


For  sufficiently  large  k,  J  f,  .=1,  hence  we  obtain  (5.4). 


(C-8) 
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